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Abstract: Dual-energy X-ray absorptiometry (DEXA) has been widely applied to measure bone mineral density 
(BMD) and soft-tissue composition of human body. However, the use of DEXA is greatly limited for low-Z materials 
such as soft tissues due to their weak absorption. While X-ray phase-contrast imaging (XPCI) shows significantly 
improved contrast in comparison with the conventional standard absorption-based X-ray imaging for soft tissues. In 
this paper, we propose a novel X-ray phase-contrast method to measure the area density of low-Z materials, including 
a single-energy method and a dual-energy method. The single-energy method is for the area density calculation of 
one low-Z material, while the dual-energy method is aiming to calculate the area densities of two low-Z materials 
simultaneously. Comparing the experimental and simulation results with the theoretic ones, the new method proves 
to have the potential to replace DEXA in area density measurement. The new method sets the prerequisites for 
future precise and low-dose area density calculation method of low-Z materials. 
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1 Introduction 

The idea of dual-energy X-ray absorptiometry can be 
traced back to the 1970s, when the method of differenti¬ 
ating low-Z materials from high-Z materials with a dual¬ 
energy scan was firstly proposed by Hounsfield [T]. Fol¬ 
lowing his work, researchers including Zatz, Alvarez et al. 
EHU reported the advantages of dual-energy method in 
material differentiation. In 1976, R. A. Rutherford uti¬ 
lized dual-energy scan to determine the effective atomic 
number and electron density of materials, and applied it 
to human tissue imaging [5]. Recently, dual-energy ab¬ 
sorptiometry has found increasingly wide applications in 
the fields of clinical diagnosis, public security, etc [6HI2]. 
As an application, dual-energy X-ray absorptiometry is 
the golden standard in bone mineral density (BMD) mea¬ 
surement, whereby percent fat in soft tissue, fat mass, 
lean tissue mass can also be obtained [MS]. 

However, low-Z materials, such as soft tissues, show 
weak absorption of X-ray. Besides, different soft tissues 
show similar characters in X-ray absorption [13] . There¬ 
fore conventional absorption-based dual-energy CT has 
an inherent defect in material differentiation for weak 


absorption materials El- While X-ray phase-contrast 
imaging offers a dramatically higher contrast for weak 
absorbers as reported in the past few decades [I5H2I], 
the reason lies in that the refraction coefficient exceeds 
the absorption coefficient with at least three orders of 
magnitude. Moreover, the refraction coefficients of dif¬ 
ferent soft tissues vary greatly m- 

In this paper, an innovative dual-energy phase- 
contrast method is proposed, which is able to remarkably 
promote the quality of identification and differentiation 
of weakly absorbing materials, and meanwhile decrease 
the radiation dose. Firstly, we derive the equations for 
area density calculation, including single-energy method 
and dual-energy method. Secondly, based on the exper¬ 
imental results and a numerical simulation, we demon¬ 
strate the feasibility of the methods. Finally, we discuss 
experimental data and potential application of the new 
dual-energy phase-contrast imaging method. 
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2 Theories 

When a beam of X-ray penetrates an object, the vari¬ 
ation of the intensity is given by 


Accordingly, we are capable of making derivation of 
the formulas to calculate the area densities with dual¬ 
energy refraction data. The phase shift of X-ray when 
passing through an object is 


I = Ioe-^^. (1) 

where Iq is the intensity of the incident X-ray, / repre¬ 
sents the intensity of the emergent X-ray, t is the trans¬ 
mission length, and /i is the linear attenuation coefficient. 
Generally, Eq. (1) is written as 

( 2 ) 

where p is the density of the substance, /i/p is known as 
mass attenuation coefficient, which is merely determined 
by the characters of the substance. And pt is the area 
density of substance, denoted by M. The area density 
thus is given by 


i^/p ■ 


(3) 


The precise measurement of BMD is crucial for the 
diagnosis of osteoporosis [22] , and dual-energy X-ray ab¬ 
sorptiometry actually measures the area density of bones 
|6|, which is primary work for further information ex¬ 
traction. Area density is usually expressed in g/cw? ^ it 
is used to describe the bone mass per unit of projected 
area, or the average mass per pixel. In order to dimin¬ 
ish the error introduced by the soft tissues, dual-energy 
method is proposed in the measurement [22l |23]. For 
each projection, two images at different levels of energy 
are obtained with the dual-energy method. For human 
body, the penetration of dual-energy X-ray can be ex¬ 
pressed by 


^ = -k J 6dt. (8) 

where k = 2'KjX is the wave vector, d is the phase factor 
and t is the transmission length. If the phase object is 
homogeneous, namely, the phase factor S is evenly dis¬ 
tributed, Eq. (8) can be rewritten as 

^ = —kSt. (9) 

Likewise, we define mass phase factor as S/p. Since we 
know that the area density M = pt^ Eq. (9) is further 
transformed to 


^ = -k{S/p)M. (10) 

As described before, the sample cylinders are even along 
length direction. Consequently, the phase shift ^ is sim¬ 
ply the function of the radius direction (we assume it as 
X direction), the relationship between phase shift ^ and 
refraction angle a is correspondingly given by 

= ( 11 ) 

Combining with the definition and assumption as is 
stated before, the area density function of the refraction 
angle is 


M = - 


[/a(a;)rfa)] 
5! p 


( 12 ) 


II — -^oLexp [—{pIp)lsMs — {p/ p)lbMb\ ■ 


Similarly, the penetration of dual-energy X-ray through 
human body can be described by 


Ih — Iqh exp [—{p/ p)hsMs — {p/ p)hbMb] • (5) 


——k{6 / p)slMs — k{6 / p)blMb. (13) 


where subscript L and H represent low and high energy, 
respectively; and subscript B and S represent bone and 
soft tissue, respectively. Supposing that Uls={p/p)ls^ 
Uhs {p/p)hsi ^lb {p/p)lb'’ ^hb {p/p)hb'i 
Rs = Uls/Uhsi Rb = Ulb/Uhs^ fhe area density of 
the soft tissue and the bone can be obtained from Eq. 
(4) and (5) as follows 


— ~k{6 / p)shMs — k{S / p)bhMb- (14) 

Supposing that = {S/p)sL^ ^sh = {S/p)sh^ ^bl = 
{^/p)bl, ^bh = {S/p)bh, Zs = ^B = " ffr ’ 
the area density of the soft tissue and the bone eventually 
become 


Mc = 


i?^ln(/^//o^)-ln(/^//o^) 
Uls ~Uhs/Rb 


Mn 


Rs / Rh ) ~ / Rl ) 

Ulb — Uhb/Rs 


(6) 

(7) 


Mb = 


^lZs-^h 
k/^BLi^B — Zs) 


(15) 


Ms = 


^lZb ~ 
k^shiZs — Zb) 


(16) 
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3 Materials and methods 


3.1 Experimental setup and working principles 


For the experimental setup, a Talbot-Lau interfer¬ 
ometer is used in combination with a tungsten rotating 
anode X-ray tube and a CCD detector constructed at 
the Institute of Multidisciplinary Research for Advanced 
Materials in the Tohoku University, Japan as illustrated 
by Fig. 1. 



Fig. 1. (color online) Layout of hard-X-ray 
Talbot-Lau interferometer, (a). The scheme of 
Talbot-Lau interferometer. (b). The working 
principle of Talbot-Lau interferometer. 


The working current of the X-ray tube is 45 mA , the 
source grating Go (period Pq = 22.7 /im, gold line height 
70 /im, size 20 x 20 mm?) is made of gold, about 80 mm 
downstream of the X-ray tube. The phase grating Gi 
(period Pi = 4.36 fim^ gold line height 2.43 /rm, size 
50 x50 mm^) is situated 106.9 mm{l) downstream of the 
source grating. The phase grating is designed to generate 
a 7r/2 phase shift at an photon energy of 27 keV. The 
samples are positioned upstream of the phase grating 
with a distance of about 50 mm. The distance between 
the analyzer grating G 2 (period P 2 = 5.4 fim^ gold line 
height 65 /im, size 50x50 mm^) and the phase grating Gi 
is 25.6 mm{d). The detector, a scintillator CCD camera 
with an effective receiving size of 68.4 x 68.4 and 

pixel size of 18 x 18 /rm^, is placed downstream of the 
analyzer grating G 2 . 

It is necessary to emphasize that the parameters of 
the optical elements are precisely designed and their rel¬ 
ative positions are well optimized so as to obtain images 
of high quality. The working principles are referred to 
[ 21 ] . The refraction angle a in the direction of axis x can 
be quantified by [24] 


A d^{x,y) 
2 tt dx 


(17) 


where x and y are the coordinates of the plane perpen¬ 
dicular to the optic axis, A is the wavelength of the inci¬ 
dent X-ray and 4>(x,^) represents the phase shift of the 
wavefront. 

3.2 Data acquisition 

Three phase objects are used for quantitative mea¬ 
surement and evaluation. One PMMA rod is 5 mm in 
diameter and two POM rods are in 5 mm and 10 mm in 
diameter, respectively. Three cylinders are all 100 mm in 
height. During the experimental operation, all of them 
are oriented parallel to the grating lines. 5-step phase 
stepping scan is adopted in the experiment. At first, the 
accelerating voltage of the tube is set at 35 kV, and the 
samples are placed at the optical axis. Then the samples 
are removed, background images at 35 kV are obtained. 
Thereafter, the accelerating voltage ramps up to 45 kV, 
and the procedure of data acquisition is exactly the same. 

3.3 Information extraction 

The absorption signal and the refraction signal are 
retrieved by the following equations [25] : 


A(m,n) = In 


k = l 
N 

E Ik{m,n) 

k = l 


(18) 


o(m,n) = X arg 


E /|(m,n) xexp(27ri|:) 


E Ik{'nT-,n) X exp{2TTi^) 


(19) 


where {m,n) denotes the pixel position of the CCD sen¬ 
sor, N is the number of the steps, subscript k represents 
that the data is obtained at k^^ step, superscript s and 
b denotes that the samples are added and removed , re¬ 
spectively. P 2 is the period of the analyzer grating G 2 , 
and d is the distance between Gi and G 2 . 

4 Results 

4.1 Information extraction 

Fig. 2 shows the absorption images and refraction 
images obtained from Eq. (18) and (19). 
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Fig. 2. (color online) Absorption and refraction 
images of the samples, (a), Absorption image at 
the tube voltage of 35 kV. (b), Refraction image 
at the tube voltage of 35 kV. (c), Absorption im¬ 
age at the tube voltage of 45 kV. (d), Refraction 
image at the tube voltage of 45 kV. The rods are 


also numbered, where is for the PMMA rod, is 
for the POM rod with the diameter of 5 mm and 
is for the POM rod with the diameter of 10 mm). 

All images are windowed for optimized appear¬ 
ance with a linear gray scale. 

Specially, Fig. 2(a) shows the absorption image at 
the accelerating voltage of 35 kV, while Fig. 2(b) is its 
counterpart of the refraction image; Fig. 2(c) shows the 
absorption image at the accelerating voltage of 45 kV, 
and Fig. 2(d) is its counterpart of the refraction image. 

4.2 Single-energy calculation 

At first, we can calculate the area density of the sam¬ 
ples with single-energy phase-contrast method. A cross 
section of the samples is selected for the calculation, 
which is denoted in Fig. 2 with red lines. The results 
are plotted in Fig. 3, and the parameters necessary for 
the calculation are listed in Table 1. 


I 








(b) (d) (f) 

Fig. 3. (color online) The results for the area density calculation along the cross section with absorption method, 
refraction method and theoretical results, which is shown with black lines, red lines and blue lines respectively, 

(a). The results of sample at the tube voltage of 35 kV. (b). The results of sample at the tube voltage of 45 kV. 

(c). The results of sample at the tube voltage of 35 kV. (d). The results of sample at the tube voltage of 45 kV. 

(e). The results of sample at the tube voltage of 35 kV. (f). The results of sample at the tube voltage of 45 kV. 


Table 1. The absorption and refraction factor at different energies and the density of the samples m- 


Tube Voltage(Mean Energy) 

PMMA 

POM 


Absorption /5(xl0 ^°) 

Refraction (5(xl0 ^) 

Absorption/3(xl0“^°) 

Refraction(5(x 10 ^) 

35 kV(25 keV) 

1.53018681 

4.26404711 

2.03051909 

5.0097696 

45 kV(27.5 keV) 

1.22495458 

3.52414105 

1.60430406 

4.1531670 

pliglcm^) 

1.19 

1.42 
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The black and red curves in Fig. 3 plot the area densi¬ 
ties calculated based on conventional absorption method 
and the proposed phase-contrast method, respectively, 
in comparison with the theoretical prediction highlighted 
by the blue curve. In order to quantify the calculation 
errors, the centric positions of the samples are chosen 
for comparison, which is also denoted in Fig. 3 with yel¬ 


low dash lines. The values and mathematical errors are 
listed in Table 2. In the table, superscript Theo and Cal 
represent theoretic values and calculation values, respec¬ 
tively, El and Eh represents the low energy (the tube 
voltage 35 kV) and the high energy (the tube voltage 45 
kV), respectively. 


Table 2. The area density comparison between the calculation values and the theoretic values at centric positions 
of the samples. 


Samples 



Absorption Method 



Refraction Method 


El 


Eh 


El 


Eh 


Kglcm?) 

%Error 


%Error 


%Error 


%Error 

PMMA 

0.595 

0.58 

-2.52 

0.57 

-4.20 

0.60 

0.84 

0.56 

-5.88 

POM 

0.71 

0.62 

-12.6 

0.69 

-2.82 

0.57 

-19.7 

0.74 

4.23 

POM 

1.42 

1.3 

-8.45 

1.19 

-16.1 

1.25 

-11.9 

1.32 

-7.04 


4.3 Dual-energy calculation 

This subsection investigates the feasibility of the 
dual-energy phase-contrast method described by Eq. 
(15) and (16) aiming to compute the area densities of 
two materials simultaneously. A numerical simulation is 
carried out on a contrast phantom whose cross section is 
shown in Fig. 4. 



Fig. 4. (color online) The cross section of the con¬ 
trast phantom used for the simulation. 


As the illustration shows, the phantom consists of a 
core cylinder made of POM, closely encircled by a cylin¬ 
drical chamber made of PMMA (r = 5 cm, R = 10 cm). 
As before, we still set the tube voltage at 35 kV and 45 
kV. Thus, the refraction angle of the X-ray 0{x) can be 
described as [23] 


0{x)= < 


‘^^PMMAX 

^/R^-x^ 


r <\x\<R 


‘^^PMMAX ‘^^PMMAX 


_^2 —x‘^ 


‘^^POmX 

—x‘2 


\x\ <r 


( 20 ) 


where we indicate the radius direction as axis x , the 
center of the cross section as origin as shown in Fig. 4. 
To simplify the simulation, we define 

0= J 0{x)dx (21) 

Fig. 5 is the schematic representation of the simu¬ 
lation process. Fig. 5(a) shows the variation of the re¬ 
fraction angle according to the pixel position at the tube 
voltage of 35 kV and Fig. 5(b) is the integration thereof. 
Fig. 5(c) and 5(d) is the corresponding simulation results 
at the tube voltage of 45 kV. 





(c) (d) 

Fig. 5. Sketch of the simulation process, (a), Sim¬ 
ulation of refraction angle at the tube voltage of 
35 kV. (b). Simulation of the integration of re¬ 
fraction angle at the tube voltage of 35 kV. (c). 
Simulation of refraction angle at the tube volt¬ 
age of 45 kV. (d), Simulation of the integration of 
refraction angle at the tube voltage of 45 kV. 
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Again the centric position of the phantom is se¬ 
lected for comparison with the theoretic values. As 
Fig. 5 shows, = 9.08765 x 10“®, 0^^ = 7.50131 x 
10“®. Based on the data in Table 1, we can figure out 
that AsL=(Vp)sL=3.58323287xl0-^ Ash={S/p)sh = 
2.96146307 x lO'^ Abl={^Ip)bl = 3.60390141 x 10“^ 
Abh={5/p)bh = 2.97178169 x 10-^ Zg = ^ = 
0.82647798, Zb = ^ = 0.82460127, consequently. 


M, 


0 r Zq — 0 f 


= 13.94(^/cm 
= 11.33b{g/cm-‘^) 


^bl{Zb — Zs) 

,, _ OlZb-Qh 

According to our assumption, however, the theoreti¬ 
cal values ought to be 


MpoM — pPOM ^ DpoM — 

MpMMA= PPMMA^ DpMMA = ^^-^{g/cm 


where Dpqm and Dpmma represents the transmis¬ 
sion length of the core cylinder and the chamber cylin¬ 
der at the centric position respectively. Naturally, the 
simulation errors are calculated as follows. 


A 

A 


]\/fCal _ J\/fTheo 

MpOM MpOM ^ = -1,83% 


l\/fTheo 
M^al _ 


IPfTheo 

^^J-PMMA 


7\TT/ieo 

XlQQ%: 


-4.75% 


where superscript Cal represents simulation values 
and superscript Theo represents the theoretic values. 


5 Discussion 


Fig. 4 shows the calculation results with the ab¬ 
sorption method and the novel phase-contrast method 
whereas Table 2 quantifies the corresponding calculation 
errors, which both proves the validity of the proposed 
method (calculation error is <5% generally). Moreover, 
the errors can result from: (1) the parameters gauge of 
the apparatus, the optical elements and the samples; (2) 
the absorption and refraction of air. Furthermore, it is 
worth noticing that there is a fearful noise jamming in the 
absorption data, conversely, the refraction data follow a 
relatively smooth curve. The distinction between the cal¬ 
culation processes is able to explain the divergence. In 
the refraction calculation method, an integration is con¬ 
tained, which is actually a filtering operation, however 
the absorption method is based on the original images, 
making the absorption method inferior to the refraction 
method in noise decrease. Nevertheless, calculation er¬ 
rors can also be introduced by the integration, further 
works is to quantify the errors. 

As stated before, the new dual-energy method is de¬ 
voted to the accurate calculation of low-Z materials, so 


two low-Z materials close to each other in effective atomic 
number {Zpmma = 6.47, ZpoM = 6.95) are employed in 
the numerical simulation. The simulation results is well 
consistent with theoretic values (simulation errors are 
—4.75% and —1.83%, respectively), and it demonstrates 
that the new method is convincing. Likewise, the error 
sources can be analyzed: (1) theoretically, the refraction 
angle at the cylinder edge approaches infinity, so the sim¬ 
ulation involves an approximation; (2) the integration to 
determine the phase shift involves errors. 

X-ray phase-contrast imaging possesses evidently 
better contrast compared to the conventional absorption- 
based imaging for low-Z materials, meanwhile they per¬ 
form similarly for high-Z materials. So the dual-energy 
phase-contrast method has certain advantages over the 
dual-energy absorption method, and they can briefly 
summarized as: (1) lower the radiation dose. With the 
novel dual-energy phase-contrast method, we can uti¬ 
lize X-ray source with lower power and lower photon en¬ 
ergy to obtain desirable images and calculation accuracy. 
Meanwhile, the radiation dose is reduced, which is criti¬ 
cal if we aim to extend its utilization to clinical applica¬ 
tion. (2) improve the calculation accuracy for bones as 
well as the soft tissues. In bone density measurement, 
soft tissue density is also measured, and the measure ac¬ 
curacy directly affects that of the bones. With the dual¬ 
energy phase-contrast method, soft tissue density can be 
measured with an improved accuracy, and it thus im¬ 
proves the bone density measurement as well. (3) reduce 
or eliminate the utilization of contrast agent in certain 
operations. In the existing dual-energy subtraction dig¬ 
ital radiography, contrast agent is requisite because of 
the similarity for soft tissues (e.g. the vessels and its 
surrounding tissues). The dual-energy phase-contrast 
method has an inherent merit in distinction of similar 
low-Z materials, so the contrast agent is not demanded 
anymore and the operation procedure is simplified. 

However, constrained by the experimental conditions, 
only two approaching tube voltages are selected (35 kV 
and 45 kV). Therefore, the influence of the tube voltage 
on the calculation precision is not fully investigated and 
understood. Besides, the phase-contrast imaging method 
we employed is based on Talbot-Lau interferometer, and 
it is generally optimized for a certain photon energy. As 
a result, the alternation of the tube voltage decreases the 
image quality. Hence we hope to develop a flexible dual¬ 
energy phase-contrast imaging method in the future. 

6 Conclusion 

In this paper, a novel phase-contrast method for area 
density calculation is presented, including single-energy 
method and dual-energy method. Experimental results 
from the refraction images and a numerical simulation 
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demonstrate that the new method is applicable to area 
density measurement. Further extension of the work in¬ 
cludes assessing the error with experiment on biological 
samples, quantifying the radiation dose, optimizing the 
dual X-ray energies and developing a flexible dual-energy 
phase-contrast imaging system. Moreover, another re¬ 
search endeavor is to apply dual-energy phase-contrast 
subtraction digital radiography in cardiovascular disease 
or inchoate cancer diagnosis, which may revolutionize the 
fields. 
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